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Abstract 

This paper is a continuation of our previous work [18]. We obtain 
two more functional relations for the eigenvalues of the transfer ma- 
trices for the sl(3) chiral Potts model at q 2 = —1. This model, up 
to a modification of boundary conditions, is equivalent to the three- 
layer three-dimensional Zamolodchikov model. From these relations 
we derive the Bethe ansatz equations. 
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1. Introduction 



One of the open problems in the theory of integrable statistical systems 
is to construct the Bethe ansatz technique for three-dimensional integrable 
lattice models. A construction of such a model is connected with the problem 
of solving tetrahedron equations [3, 4] which insure integrability of a three- 
dimensional model. These are a system of thousands of equations in the 
simplest nontrivial case. Hence, the problem of solving them is very difficult. 

There are only a few known integrable three-dimensional models which 
are interesting from the physical point of view. The first nontrivial exam- 
ple of such a model was proposed by Zamolodchikov in 1980 in [1,2]. The 
tetrahedron equations for the Zamolodchikov model were proved by Baxter 
in [5]. 

Bazhanov and Stroganov [6] observed that the Zamolodchikov model and 
the three-dimensional free-fermion model were "weakly equivalent", i.e., the 
free energy of the Zamolodchikov model and the free-fermion model satisfied 
the same symmetry and inversion relations. The assumption that analytical 
properties of the free energy were also the same resulted in a coincidence 
of the free energy for the Zamolodchikov model and the free-fermion model. 
In 1986 Baxter [7] calculated the partition function for the Zamolodchikov 
model with some modification of boundary conditions for the case of the 
infinite cubic lattice and for the case of the lattice which is infinite in two 
dimensions and finite in the third one. His result was similar to the result 
by Bazhanov and Stroganov for the free-fermion model but not the same. 
Namely, the partition function for the Zamolodchikov model was made up 
of a sum of two parts. The first part coincided with the partition function 
of the free fermion model and had the usual analytical properties for two- 
dimensional models. The second part was expressed in terms of the Euler 
dilogarithm function and had the cut in the complex plane. Therefore, the 
assumption that the free energy for the Zamolodchikov model and the free- 
fermion model had the same analytical properties was incorrect. However, 
the similarity of these results was remarkable. Later Baxter and Quispel in 
[17] tried to clarify this fact. Namely, they constructed the Hamiltonian for 
the two- and three-layer Zamolodchikov model. The two-layer case turned 
out to correspond to the two-dimensional free-fermion model. The Hamilto- 
nian for the three-layer case contained cubic interaction terms and seemed 
not to be the Hamiltonian for the free-fermion model. 
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Another important step in the theory of the integrable three-dimensional 
models was done by Baxter and Bazhanov in 1992. Namely, they observed 
[9] that the sl(n) chiral Potts model at q 2N = 1 [14, 15] was equivalent to the 
n-layer three-dimensional model which turned out to be the iV-state gener- 
alization of the Zamolodchikov model. It was also mentioned that as for the 
Zamolodchikov model this equivalence is valid only up to some modification 
of boundary conditions which should not effect the partition function in the 
thermodynamic limit. The partition function for the Baxter-Bazhanov model 
was calculated in their next paper [10]. The result appeared to be connected 
in a remarkably simple way with that for the Zamolodchikov model. 

We hope that a development of the Bethe ansatz technique for the Zamo- 
lodchikov and Baxter-Bazhanov models could shed a new light on the prob- 
lems discussed above. Since the n-layer case of the Zamolodchikov model is 
equivalent to the sl(n) chiral Potts model at q 2 = —1, up to a modification 
of boundary conditions, we can try to construct a Bethe ansatz for the sl(ri) 
chiral Potts model. 

The Bethe ansatz technique is usually applicable to the study of effects 
connected with the finite size of a lattice. Therefore, there is a good chance 
that it will be useful for an investigation of the finite size corrections and the 
excitations. 

Our first step is to develop this program for the three-layer case of the 
Zamolodchikov model with modified boundary conditions, i.e., for the s/(3) 
chiral Potts model at q 2 = — 1. 

This work is a continuation of our previous paper [18] where some func- 
tional relations for the si (3) chiral Potts model at q 2 = —1 were derived 
and the nested Bethe ansatz was constructed in the particular case when the 
vertical rapidity parameters coincide. Unfortunately, we did not succeed in 
solving these functional relations. Our goal here is to derive other functional 
relations and to obtain from them the Bethe ansatz equations for the general 
case. 

The paper is organized as follows. In Section 2 we recall the basic formu- 
lations of the sl(3) chiral Potts model and its correspondence to the modified 
three-layer Zamolodchikov model. In Section 3 we fix the definitions of the 
transfer matrices and discuss some of their simple properties. In Section 4 
we give two functional relations for the eigenvalues of the transfer matrices. 
In Section 5 we obtain the Bethe ansatz equations. In the last section we 
give a brief discussion of the obtained results and directions for further in- 
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vestigation. In the Appendix we outline the basic steps of the proof of one 
of the functional relations. 



2. Basic formulations. 



The basic formulation of the Zamolodchikov model and it's generalization, 
the Baxter-Bazhanov model can be found in papers [1,2] and [9]. In the 
last paper it was observed that the Boltzmann weights for the sl(n) chiral 
Potts model at q 2N = 1 were a product of the n more simple weights (see 
formulae ( 3.7 - 3.13 ) of [9]). Hence, the "star" weight for sl(n) chiral Potts 
model appeared to be a product of the n weight functions interpreted as the 
Boltzmann weights for some TV-state three-dimensional model: 



(6i...6 n ) 





Figure 1. 



Each weight of the product in the RHS depends on the eight spins with 
iV possible values. For the case N = 2 this model turned out to be just the 
Zamolodchikov model rewritten by Baxter in the "interaction-round-cube" 
form [5]. As it was already mentioned in Introduction this equivalence is 
valid up to some modification of boundary conditions. 

Since we study the three-layer case of the Zamolodchikov model we need 
to consider the s/(3) chiral Potts model at q 2 = — 1. The basic notations 
of this model were adduced in papers [14, 16] or in [9]. It can also be found in 
our previous paper [18] but to be independent here we give some necessary 
basic definitions below. 

The model is formulated on the square lattice (see Figure 2) 
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The interaction is defined by two types of the weight functions W pq (a, (3) 
and (W qp (a, P))" 1 which depend on the neighbouring spin variables and spec- 
tral parameters. The rule how to choose these weights is shown on Figure 3: 



V 



p / 

/ 



p 
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W qp (a,(3) 
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Figure 3. 



The Boltzmann weights depend on the rapidity parameters. Each rapidity 
variable is represented by three 2- vectors {hf(p), h~(p)), i = 1,2,3 which 
specify the point p of the algebraic curve T defined by relations 



K(pY 



Wi,j = 1,2,3, 



(2.1) 
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where Kij are 2x2 complex matrices of moduli satisfying 

detKij = 1, K« = KijKjkKkt = 1 (2.2) 

and indices k take values 1, 2, 3 modulo 3. 

Further we will need the automorphism r on the curve T defined as follows 

h+(T(p)) = h+(p), hj(r(p)) = -hj(p), j = 1,2,3. (2.3) 

The curve T can be defined in a different way which is also useful. Namely, 
for two arbitrary points p and q on this curve the following combination 

Apq = ht{pfh-{ q f - h-(p) 2 ht(qf, (2.4) 

should be the same for all i = 1,2,3. It is easy to see that both these 
definitions are equivalent to each other. 

The Boltzmann weights depend also on spin variables. Each spin variable 
is described by a two- vector 

a = (ai, a 2 ), oii G Z2 2 = 1,2. (2.5) 

Then the function W pq (a, (3), a, (3 G Z 2 x Z 2 is defined as 

W p , q (a, (3) = (-l) Q Mg„(0, a - (3), (2.6) 

where 

Q{a, (3) = A (A - ai ) + fo(Pi -a,+p 2 - a 2 ), a,(3eZ 2 xZ 2 (2.7) 
and the function ^(0, a) has the following form 

"it \hf(p)h;( q )-hf( q )h;(p)(-l)-P) 
9 Pq (0, a) = . (2.8) 

n n (hnp)K( q ) - 

We choose a normalisation of W pq (a,(3) as 

W pq (0,0) = l. (2.9) 
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Then it is easy to see that 

W pp (a,(3) = 5 a ^ (2.10) 

where 

g a = P (mod 2); 

01,13 1 0, otherwise. 

The function W pq (a,{3) satisfies the inversion relation 

]T W pq (a,P)W qp (P n ) = 5 a ^ pq , (2.12) 

where the inversion factor $ pg is given by: 



% q = ^ (2-13) 



and 

3 3 

^ = n^(P)^"(9) + UK(p)K(q), (2.14) 



i=l i=l 
3 



D Pq = U(ht(p)K(q) + h-(p)ht{q)). (2.15) 
i=i 

As it was shown in [8] (see [9] for details) the Boltzmann weights W 
satisfy the "star-star" relation which provide the integrability of the sl(n) 
chiral Potts model. This relation looks as 

Wp> p (6,a) . Wg^ptJ) 

W plp (j, (3) qq 1 ' 1 1 7 ' ] qq 1 ' 1 1 7 ' } W q , q {5, 7 ) ' 1 } 

where two "star" weights are defined as follows 

<K ft 7, <5) = E Sii^S^ffi^ , (2 . 17) 



TF p ^(a,/5) 

The objects defined in (2.17-2.18) are the "star" -weights. As it was men- 
tioned above these weights correspond to the three-layer case of the Zamolod- 
chikov model. To be exact for the general case of the rapidity variables hf(p) 
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satisfying (2.1) the corresponding Zamolodchikov model is inhomogeneous in 
the third direction. In fact, we will be interested mainly in the homogeneous 
case: 

ht(p) = 1, K(p) = p. (2.19) 

It is easy to see that the defining relations (2.4) are trivially satisfied. There- 
fore we do not need to work with the high genus curve T. 

As it was pointed out in [10] the rapidity variables can be parameterized 
in terms of the spherical angles and excesses #i, #2, 



1 



—i tan 



0, 



p_ 
p> 



% tan 



p 

q 



e-^Wtan-tan-. (2.20) 
V 2 2 



3. Transfer matrices 

Here we use slightly different definitions of the transfer matrices comparing 
with [18]: 



k=i 



W q/q (i k ,i k+1 ), 



Tin- a n '\ jl '"^ N - TT W q'qUkJk+l) W pq/ (j k+1 , i k ) 

1 \P) <j, y )i lt ...,i N 11 



k=i 



W pq (j k ,i k ) 



which are shown on Figures 4 and 5: 



(3.1) 
(3.2) 



T(p; q,q>) 




ll / \ 12 / \ «3 ■■• «JV / \ ll 

q' q q' q q' q 

Figure 4. 
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h \ / h \ S ja ■■■ jN \ S ji 



T{p; q, q') 




Figure 5. 

where we imply the cyclic boundary conditions in+i = H an d Jjv+i = ji- 

We note that these definitions differ from the previous ones just by the di- 
agonal equivalence transformation. Of course, it does not effect the partition 
function. 

Below we shall use more simple notations T p = T(p; q, q') and T p = 
Tip; q, q') assuming that the rapidities q and q' are fixed. Due to (2.16) these 
transfer matrices T p and T p commute. Namely, for two arbitrary rapidities 
p and p' 

[T p ,T p/ ] = [T P ,T P ,] = [T P ,T P ,] = 0. (3.3) 

One can consider some limiting cases. From (2.6-2.9) we can conclude 
that if q' ^ p we have 

T p = X~\ T p = X, (3.4) 
where X is the shift-operator: 

N 

k=i 

if q — > p then 

T p = I, (3.6) 
while T has the singular matrix elements. 

4. Functional relations 

Further we will consider only the case of the homogeneous three-layer Zamo- 
lodchikov model. Due to the commutation relations (3.3) we can diagonalize 
the transfer matrices T p and T p simultaneously. 
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Let us denote eigenvalues of T p and T p by t{p) and t(p) where we omit a 
dependence on q and q'. 

In Appendix we outline the proof of the following couple of functional 
relations 

t(p) t(p) t(-p) t(-ujp) = 
<j%t{p)t(-up) + <fil(-p)t{-up) + <$t{p)t(up) + $ l(-p)t{up) 

(4.1) 

and 

t(-cop) t(-p) t(p) t(p) = 
(f>'o t(-up)t(p) + <j>'?t(-up)t(-p) + 0f t{up)t{p) + 0f t(cup)t(-p), 

(4.2) 

where 

(p + wg) (p + c^g) (p + wg') (p + w - Y) 

00 = 4 bT^p ' 01 = 4 HTW ' (43) 

, , (P - q) (P + ^qf (P + uj~ l q') (p - uq') 
<p2 = 4 



(p — u> 2 g) (p + q) 2 (p + g') (p — uj 1 q') 
= 1 (p - q') (p + u'V) 2 (p + uj~ l q) (p - ^g) 
(p - cj 2 g') (p + g') 2 (p + g) (p - w _1 g) 

and 0- can be obtained from 0j by the substitution g — > — g. Here uj is the 
root of unity of power three 

uj = e s . 

From the limiting cases (3.4) and (3.6) we have some initial data: 

t(p;q,p) = n, t(p;q,p) = SI' 1 , t(p;p,q') = l (4.5) 
where fl is some root of unity of power N: 

n N = i. (4.6) 

From (4.1) and (4.2) one can see that the pair of functions t' and t: 

t'(p; q, q') = t(p; -q, q'), t'(p; q, q') = t(p; -q, q') (4.7) 
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satisfy the same relations (4.1) and (4.2). However, it is not true that 
t(p; q, q') = t(p; —q, q') for all eigenvalues. The transformation (4.7) inter- 
changes also the eigenvectors of the transfer matrices which belong to the 
same symmetry sector. 

The analysis of the eigenvalues t(p) and t(p) shows that it is convenient 
to extract some "kinematic" multipliers: 

2 n _ 2 N 

= ¥TW¥TW S(P) - t(r) = ¥^W¥TW" (P) (4 ' 8) 

where s(p) and s(p) are the polynomials of the degree n in the variable p. So 
far, we have no proof that the degrees of s(p) and s(p) should be the same. 
Therefore, we accept it conjecture. 

Substituting the definitions (4.8) into (4.1-4.2) we obtain the functional 
relations for s(p) and s(p): 

s(p) s(p) s(—p) s(—up) = 
K «(p) s(-ujp) + Af s(-p) s(-ujp) + A^ s(p) s(ujp) + A^ s(-p) s{up) 

(4.9) 

and 

s( — up) s(—p) s(p) s(p) = 
X'o s(-up) s(p) + A'f s(-up) s(-p) + X'z s(ujp) s(p) + A'f s(up) s(-p), 

(4.10) 

where 

A = (p +uq) (p +W 1 q) (p + q') (p - q'), 
\ l = (p +uq') (p +uj- 1 q') (p + q) (p - q), 
A 2 = (p - q) (P +^~ 1 q) (P - uq') (p - q'), 
A3 = (P ~ q') (P +uj' 1 q') (p - uq) (p - q) (4.11) 

and A- can be obtained from \ by the substitution q — > —q. 
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tip) 



5. Bethe ansatz equations 



To construct the Bethe ansatz we consider zeros of the polynomials s(p) and 
s(p): 

n n 

s(p) = a n (q,q') J[ (p - p t ), s(p) = a n (q,q') ]J (p - p { ) (5.1) 

where the power n takes only two possible values 2N and 2iV — 1. The 
functions a n and a n should be compatible with the initial conditions (4.5). 
Unfortunately, it is not easy to calculate them explicitly but their product 
looks very simple: 

a 2N (q, q') a 2N (q, q') = 4, a 2N -i{q, q') a 2 N-i(q, <f) = N (q' 2 - q 2 ). (5.2) 

Now we can set p to be some zero of the LHS of (4.9) and consider the 
equations which follow from the RHS. In fact, we have four possibilities to 
do this: 

P^Pi, p->-Pi, p-t-w^Pi, P^Pi- (5.3) 

It is not difficult to obtain that first three possibilities give us two different 
sets of Bethe ansatz equations: 



f(Pi,u,-q') j= iPi-^ + Pj 



and 



where 



fjp^, g) jy = ( -l ) n-lfl ft+^ (5 . 5) 

f(p,x,q) = - — (5.6) 
p + q 

The fourth possibility in (5.3) gives some complicated compatibility condi- 
tions for the solution to the Bethe ansatz equations (5.4-5.5). Of course, Pi 
and p { are the functions of q and q' . A similar consideration of the second 
functional relation leads to the same Bethe ansatz equations (5.4) and (5.5). 
It is obvious that s{p) and s{p) are homogeneous in p, q, q' . So let 

q = 1, p = ix, q' = iy (5.7) 
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where x, y are real. 
Conjecture 

s(x,y) = s*(x,y) (5.8) 
We checked it numerically for N = 2, 3. Let us set 

Pi = iri{y), Pi = i7i(y) (5.9) 

Then 

n n 

s(p) = a{y)i n ]J(x- n(y)), s(p) = a{y)i n ]J(x - n{y)) (5.10) 
i=i i=i 

From (5.4) we obtain 

a(y) = (-l) n a*(y), r l (y)=r*(y), i = l,...,n (5.11) 

It is easy to see that (5.5) can be obtained from (5.4) by a complex con- 
jugation. Somehow this is the "proof" of the conjecture (5.8). Then we 
have 

^-#'Mi» =( . r ^.+^ i e = ±1 . (5 . 12) 



(n + uj £ y) {in - 1) 



j=1 n-u *r* 



One can obtain from (5.12) the set of equations on absolute values and phases 
of n. 



6. Discussion 

In this paper we have only presented the Bethe ansatz equations. We shall 
give the detailed analysis of these equations elsewhere. The technique we 
use here is in the spirit of the Baxter Q-matrix method [12]. The role of the 
Q-matrices is played by the one-layer transfer matrices. It corresponds to the 
result obtained by Bazhanov and Stroganov in [11] for the chiral Potts model. 
We think that the algebraic Bethe ansatz technique can also be developed. 
However, there are some problems like an appropriate choice of the reference 
state which are out of our understanding so far. 

We should note that the functional relations we have derived here and 
those which were obtained in [18] can be considered together. Perhaps the 
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combining of all these relations could give more information about the eigen- 
values t(p) and tip). 

We hope that the result obtained by Baxter in [7] for the partition func- 
tion of the Zamolodchikov model on the lattice oo x oo x 3 can be reproduced 
in the thermodynamic limit of the Bethe ansatz equations (5.4-5.5). We also 
hope that a standard program of a study of the excitations and finite size 
corrections 3 can be performed. 

We think that the technique described in Appendix can be generalized 
to the sl{n) case. In principle, a general procedure seems to be more or less 
clear. However, the technical difficulties one could face can be, of course, 
much more serious. 
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Appendix 

Here we outline the derivation of (4.1). The second relation (4.2) can be 
obtained in a similar way. In fact, we derive it for a general case of the inho- 
mogeneous Zamolodchikov model. Our key relation for the transfer matrices 
T p and T p we would like to obtain looks as follows: 

T p T p T T ( P ) T p * = &q Tp T p * + <3>f T r(p ) T p * + $^ T p T r(p *) + $^ T r(p ) T r ( p *) 

(A.l) 

where p* is one of two nontrivial solutions of the equation 

H t Hi /A , 

7F = "TP (A ' 2) 

Hp rip* 



3 See for example the book [13] and references therein. 
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where 



3 AD PA 

and 

$ = 4-^-, $! = 4-^1, $ 2 = 4« M v, $3 = 4 u M Up,/. (A.4) 

The function i?,/} and A are given by the formulae (2.14,2.15) and (2.4) 
respectively. 

In fact, in our previous paper [18] we made the first step. Namely, we 
expressed the matrix product T p T T ^ as a sum of two terms 4 . The first 
one corresponds to the first term in the RHS of (A.l). The second term 
was written in terms of some L-operators. When the vertical rapidities q 
and q' coincide this L-operator corresponds to the second fundamental rep- 
resentation 3 of the quantum s/(3) algebra. Therefore we shall denote it as 
L. 

Now we have to do the next step. Namely, we should consider the matrix 
product: 

,T T J"! 

i=i 



N 



(T P L P )\;l = Trl[B^ i+1 (q',q;p), (A.5) 



where 



L is given by 



\BIM^V\, - a), (A.6) 



W (B m) 

LiA^a) = £ n) W T{p)q ,(a, n) W q/p (n, (5) pq{fJ / ) C(j, m), 

(A.7) 

and all indices a, (3, , ~f,5,n,m are two-component vectors taking one of four 
possible states (0, 0), (0, 1), (1, 0), (1, 1), i, j = 1,2,3, 

4 Wc considered two cases A = 0, 1 corresponding to two automorhisms t a (see formula 
(3.2) of [18]). Here we consider only the case A = 1. A consideration of the case A = 
gives the same result. 
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C(2h + k 2 ,n) = l(_ 1 )fcmi + fc 2 ™ 2 _ ( A8 ) 

Inserting the identity matrices 3x3 between each pair of B in the RHS 
of(A.5) 

3 

1 = x <M^ a )' ( A - 9 ) 



i=l 



where 



L (l,a) = (1,0,0), 

L (2,a) = (-(-ir +a2 ,l,O), 

L (3,a) = (-(-ir,O,l), (A.10) 



i? (i, a ) = (i,(-ir +a2 ,(-ir; 

fl (2,a) = (0,1,0), 

<M3,a) = (0,0,1), (A.ll) 
one can check that the transformed matrices: 

[BU^^a = Mhl)Bl 5 (q',q;p)M3,5) (A.12) 

satisfy the following property 

fc(?',?;p)] 21 = K 5 (g',g;p)] 31 = (A.13) 

for all possible values of indices a, 7, 5. 

Therefore, we have a decomposition 1 + 2. It is not difficult to check that 

fete',,;^, = V (-ir**^g$ (A,4) 

where $ pg is defined in (2.13). In the RHS of formula (A. 14) we can recognize 
the "building block" of the transfer matrix T r ( p ). So, after taking the product 
and trace as in the RHS of (A. 5) we obtain the second term in (A.l). 
Now let us define 2x2 matrices with elements 

K s (q',q;p)]. 3 = K s (q',q;p)]. +1 . +1 , i,j = 1,2. (A.15) 
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The matrices B^ 5 (q', q;p) have the form of the following matrix product: 
B° s tf,q',p) = V^(a, 7 )l/ M (M) (A.16) 

where 

[U Pq (a,S)] tj = ExL(i;m,a) t^w^ ( M^n,5) (A.17) 

n,m ''PI 

[171,0 )W T (p) q {a,n) 

and 

[Vpq(a, 5)] . . = Xl(i; m, 5)W r{p)q (a, m)W qp (m, n)W pq (n, 5)xr{j] n, a). 



n,m 



(A.18) 

Here we use the following notations 

3 

XHR)(i;m,a) = ^C{k,m) [0 L (H)(i + l;a)] fc , i = 1,2. (A.19) 
fc=i 

It is interesting to note that U pq and V pq satisfy the property which is 
similar to that for B given by (A. 13) 

[U Pq (a,5)] iQ = [V pq (a,5)} t0 = 0. (A.20) 
In addition, we have 

[U Pg (a,5)] 00 = -^(-ir + ^W qr(p) (5,a) (A.21) 



[V Pq (a, S)] m = (-ir +h W T{p)q (a, 5) (A.22) 



~4 

and $ pg is given by (2.13). 

Using the definitions (A.17) and (A.18) we obtain 



= (-l) ai+aa 77 M (ai,a 2 ;(Ji,(J 2 ) 
=(-l) ai+5 ^ p(? (a 1 ,a 2 ; ( 5 1 , ( 5 2 ) I z ^ q . a ^ 5) 



Z32(p,q;a,5) 1 /' 

(A.23) 



2 *7p,(ai,a 2 ;5i,5 2 ) ( _ ^ i_ 

72 (p,<2) 



(A.24) 
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where 

^(p.j.a.i) = (-1) Mptq) _ ( _ l)n+Sl)rMq) • (A.27) 

*,<p, = (-1) (73(M) _ (_ 1 )o,+Ji+o,+&) 72 (p i9 )- < A ' 28 > 

It is easy to see from (A. 23) and (A. 24) that the matrix U pq is connected 
with V pq by the matrix inversion up to some coefficient: 



V pq (a, P)U pq (a, (3) = (71 (p, q)-f 2 (p, q)lz{p, q) - 1) x 
(72b, g) - (-l) Q2+52 )77p, ? («i, a 2 ; 5 1 , 5 2 ) 2 



(A.29) 



(7i(P,9) - (~l) ai+Sl ) (7s(p,5) - (-l) ai+5l+a2+<52 ) 7 2 (p,g) 2 
The important fact is a degeneration of these matrices which occurs when 

7i(P» 9) 72(P, 5) 7s(p, 5) = !■ (A.30) 

Using the "star-star" relation for the Boltzmann weights and property 
(A. 20) we can prove that the matrices U and V should satisfy the following 
important relation: 

]T V^(a, 7 )C^(«> S ) V„,(6',6) W^tf, a)W p , q (a, 5') = 



E W p/q ( r p)W q , p ,((3,6) V pp ,{i, 7)^(7', P)V P A6',(3). 

(A.31) 

From this relation we can deduce that choosing the rapidity variable p' 
in such a way that all matrices V pp i are degenerate we get the decomposition 
of the matrices with the following elements: 
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Wf(q',q;p,p')] tj = J2ly P A^l)U pq (a,5)} ij W q>p ,{i, a)W p , q (a, 5'). 

(A.32) 

It means that we can reduce the matrices Dj^f (q',q;p,p') by the quasi- 
equivalence transformation to the upper-triangular form. This technique is 
rather similar to that which was used by Baxter for a derivation of the Q- 
matrix equation for the 6- vertex and 8- vertex models [12]. 

So, first we should choose the point on the curve p' to provide the degen- 
eration of matrices V p y. Therefore, we should fulfil the condition (A. 30) for 
the pair (p,p'): 

Ii(p,p')l2(p,p')l3(p,p') = 1- (A.33) 

This equation has three solutions (up to some choice of signs). One of 
them 

A pp , = (A.34) 
corresponds to the automorphism r: 

p' = r(p). (A.35) 

Two another solutions can be obtained by taking the second power of (A.33) 
and using (2.1). In this way we arrive to the quadratic equation for the 
coordinates of p' with coefficients depending on coordinates of p. Let us 
denote its roots as 

P± = r±(p). (A.36) 
Let us choose one of these solutions, for example, 



P* = P+ (A.37) 

and set the point p' in the formulae above to be p*. 

It is easy to conclude from (A. 24) that up to some factor the matrices 
V pp * are proportional to 



Vpp*(a, 5) 



( 



(-1)«2( 7 *_(--1)"1+< 5 1) 



1 V»a ^2-(- l ) ai+Sl+a2+S2 



"7| 



(A.38) 
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where 



7* = liipy)- (A.39) 
So, the vectors which provide the decomposition 1 + 1 can be chosen as: 



Cl(1;M) = (1,0), ( L (2;a,5) = (-(-1)°* 7l72 ^ , 1), 



(A.40) 



U2;a,5) = (0,l), C R (l;«,5) = (l,(-ir^ r V 1 ' ' 



7l72 ~~ ( — l) a i+ 5 i+ a 2+<5 2 y ' 

(A.41) 

which satisfy the natural condition 

X)Cfl(i;M)xa(i;M) = (A.42) 

Now let us consider the transformed matrices IX: i : 

K/V,9;p)] y = Ui\i,l)^fti,<l\™*)tej\t',S) (A-43) 

where i, j = 1, 2. 

One can check the decomposition property: 

[5?/ P)] 21 = 0. (A.44) 

Now we should study diagonal elements of these matrices [D^f (q 1 , q; p)\ .., 
i = 1,2. 

It can be checked that the following expressions for Da are valid 

[^/V,?;P)]« = A ^ (^^^P)|^^ ) (7,5,7 / ,5 , ), (A.45) 

where W ? (? p / (? t{p * ) (7, 5, V, 5') is given by (2.18) and 

Pi = P, P2 = r(p). (A.46) 
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For the scalar functions A, we have 



Ai(g',g;p) = 4:V q ,, p u q , p , A 2 (g',g;p) = 4:U q/ , p v q , p 



(A.47) 



and the functions u and t> are given by (A. 3). The gauge matrices Ai are 
given by: 



A x = 



I 


1 


ai 


ai 


1 


\ 




d\ 


1 


1 


ai 








-1 


-1 


— ai 




\ 


-1 


— Oi 


— ai 


-1 


/ 


I 


1 


a 2 


— a 2 


-1 


\ 




-a 2 


-1 


1 


a 2 






a 2 


1 


-1 


— a 2 




\ 


-1 


-a 2 


a 2 


1 


/ 



(A.48) 



(A.49) 



where 



ai 



^(P)^(P*) + ^3 + (p)^3"(^) 
Kip) h 3 (p*) - ^(p)^(p*)' 



a 2 



-1/ai. 



(A.50) 



So, we have succeeded in reducing the four-index objects [D^f iq', Q',p)] u 
to the original "star"-form. It is not difficult to observe that after taking the 
product and trace we obtain the last two terms in (A.l). Using commutation 
relations (3.3) we can simultaneously diagonalize the transfer matrices and 
get the functional relation for the eigenvalues. Taking into account that for 
the homogeneous case of the Zamolodchikov model the automorphism r acts 
just as negating of p: 

rip) = -p. (A.51) 



and rapidity p* can be taken as 



p = —cop 



(A.52) 



we come to (4.1). 
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